Accounting for Baryons in Cosmological Constraints from Cosmic Shear 
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One of the most pernicious theoretical systematics facing upcoming gravitational lensing surveys 
is the uncertainty introduced by the effects of baryons on the power spectrum of the convergence 
field. One method that has been proposed to account for these effects is to allow several additional 
parameters (that characterize dark matter halos) to vary and to fit lensing data to these halo 
parameters concurrently with the standard set of cosmological parameters. We test this method. 
In particular, we use this technique to model convergence power spectrum predictions from a set 
of cosmological simulations. We estimate biases in dark energy equation of state parameters that 
would be incurred if one were to fit the spectra predicted by the simulations either with no model 
for baryons, or with the proposed method. We show that neglecting baryonic effect leads to biases 
in dark energy parameters that are several times the statistical errors for a survey like the Dark 
Energy Survey. The proposed method to correct for baryonic effects renders the residual biases in 
dark energy equation of state parameters smaller than the statistical errors. These results suggest 
that this mitigation method may be applied to analyze convergence spectra from a survey like the 
Dark Energy Survey. For significantly larger surveys, such as will be carried out by the Large 
Synoptic Survey Telescope, the biases introduced by baryonic effects are much more significant. 
We show that this mitigation technique significantly reduces the biases for such larger surveys, but 
that a more effective mitigation strategy will need to be developed in order ensure that the residual 
biases in these surveys fall below the statistical errors. 

PACS numbers: 98.80.-k,98.62.Py,98.35.Gi 
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I. INTRODUCTION 



Weak gravitational lensing is a potentially powerful 
probe of cosmology (e.g., Refs. [ll-Q) ^- Imaging sur- 
veys such as the Dark Energy Survey (DES) and, in the 
longer term, the surveys of the Large Synoptic Survey 
Telescope (LSST), the European Space Agency's Euchd 
satellite, and the Wide Field Infra-Red Survey Telescope 
(WFIRST) expect to measure the power spectrum of cos- 
mological weak lensing with sufficient precision to im- 
prove constraints on dark energy dramatically. However, 
a number of sources of systematic error must be con- 
trolled in order to achieve these goals. From a theoretical 
perspective, it is necessary to predict matter power spec- 
tra with precisions of better than one percent over a wide 
range of scales [§, [l3|- This is a challenging goal, but 
significant progress has been realized utilizing N-body 
simulations containing only dark matter [ill- 14 1. The 



largest remaining challenge to this goal is to account for 
the influence of the baryonic component of the universe 
in these predictions. Baryonic effects have been shown to 
alter lensing power spectra significantly on small scales 
[15I - I20I . This theoretical systematic error associated with 
baryonic processes is sufficient to cause large systematic 
errors in inferred dark energy parameters if unaccounted 
for [l3,[23-il, though Ref. UM explored methods to cull 
data in order to protect against scale-dependent uncer- 
tainties in predicted power spectra. In the present work, 
we assess a proposal to mitigate dark energy biases in- 
duced by baryonic efi^ects using a method proposed in 
Ref. 



21 



^ This application of lensing goes back more than forty years (e.g., 
Ref. [1) 



Rudd et al. [18| recognized changes in the internal 
structures of dark matter halos as the cause of the largest 
alterations to lensing spectra in baryonic simulations (a 
result confirmed in Refs. [l^ H^)- Consequently, Zent- 
ner et al. HJ] suggested a strategy to mitigate baryonic 
effects in forthcoming lensing analyses. Zentner et al. 
[2l| proposed altering the canonical relationship between 
halo mass and halo concentration (e.g., [U, [25|) to ac- 
count for the matter redistribution precipitated by bary- 
onic effects, as this enables the simulations of Ref. UM 
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to be modeled successfully. Zentner et al. [2l[ then sug- 
gested introducing additional parameter freedom into the 
concentration-mass relation, a necessity because this re- 
lation cannot be unambiguously predicted due to the un- 
certainties in baryonic processes, and fitting the data si- 
multaneously for the parameters that quantify the mass- 
concentration relation and the cosmological parameters. 
The value of this strategy is that it can reduce systematic 
errors (or biases) in dark energy parameters to acceptable 
levels, while increasing the statistical errors on dark en- 
ergy parameters by only ~ 10 — 40%, depending upon the 
experiment and the complexity of the mass-concentration 
relation plj (a similar argument can be made for modi- 
fied gravity [23]). 
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Parameter 1 

FIG. 1: Cartoon view of the effect of bias on cosmological 
parameters (here called "Parameter 1" and "Parameter 2"). 
The true value of the parameters is given by the red star. If 
the bias (in our case the effect of baryons on the weak lensing 
power spectrum) is not accounted for, the allowed region in 
parameter space will be given by the shaded blue region at the 
top right. The parameters will be offset from their true val- 
ues, or biased. We call the offset in this case the "raw bias." If 
one attempts to mitigate the bias by introducing new param- 
eters (in our case allowing for a varying mass-concentration 
relation), the allowed region will shift to that given by the 
shaded green contour at the lower left. The errors are now 
larger due to the increased number of parameters used in the 
fit, but the offset, which we refer to as the "residual bias," is 
much smaller than the raw bias. 

Our aim here is to test this mitigation strategy more 
extensively. We wish to determine if this algorithm will 
extract cosmological parameters successfully from up- 
coming survey data. Successfully here has a specific and 
technical meaning, a cartoon version of which is illus- 
trated in Fig. [1] First, success demands that biases in 



the cosmological parameters due to inaccuracies in the- 
oretical models should be small. There are two biases 
at play here, the "raw bias" (the offset of the smaller 
contour in Fig. [1} before any mitigation is applied and 
the "residual bias" (offset of the larger contour in Fig. [T]) 
which remains after fitting for the new free parameters. 
Ideally the residual bias will be much smaller than the 
raw bias; for the method to be truly effective, the residual 
bias should be smaller than the statistical error. Second, 
success requires that the additional parameter freedom 
introduced by the model should not inflate the error bars 
on cosmological parameters so much as to markedly re- 
duce the constraining power of the experiment. At min- 
imum, the increase in the statistical error bars due to 
additional parameters should not be so large as to nullify 
the reduction in the systematic errors. 

To carry out this test, we use the results from the Over- 
whelmingly Large Simulations (OWLS) 0, M, ^ as 
mock data. The OWLS suite consists, in part, of a set 
of ten simulations, each with the same initial conditions 
evolved in the context of the same cosmology. One simu- 
lation treats only dark matter, while the other nine model 
baryonic processes using different effective models. We 
proceed by assuming that each one of the OWLS sim- 
ulations, in turn, produces the true matter power spec- 
trum. We fit each of the OWLS predictions for lensing 
power spectra with our mitigation model, including nui- 
sance parameters. We compute residual differences in 
power spectra between our fits and the OWLS predic- 
tions and use these differences to estimate the biases in 
dark energy equation of state parameters that would be 
realized after applying the mitigation scheme. We re- 
peat this analysis in the context of two distinct imaging 
surveys. The first survey we consider has the precision 
expected from DES. The second survey we consider rep- 
resents more long-term. Stage IV^ surveys, such as may 
be conducted by LSST or Euclid. 

We will show that baryonic effects may reasonably lead 
to raw biases as large as ^ 2cr — 6cr (where a represents 
the marginalized statistical error) on dark energy equa- 
tion of state parameters if unaccounted for in the analysis 
of DES-like data. The size of the bias depends upon the 
range of multipoles used in the analysis and the baryon 
model. This broadly confirms prior estimates |20l - [22 |. 
We will then show that this mitigation scheme can render 
systematic errors sufficiently small, so as to suggest con- 
centration fitting as an attractive strategy for the cosmo- 
logical analysis of lensing power spectra from DES. In all 
cases that we consider, the residual biases remain < 0.5f7 
and can be kept < 0.1 cr if the range of scales included 
in the cosmological analysis is restricted to ^ < 2000, 
though restricting scales comes at a non-negligible cost 
in statistical error. 



^ Using the classification scheme of the Dark Energy Task Force 
[Ssll . within which DES would be a Stage III experiment. 
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For Stage IV experiments with wide sky coverage, such 
as LSST or Euchd, the conclusion is shghtly more com- 
pUcated. Absent any mitigation scheme for baryonic pro- 
cess such future experiments may be subject to raw biases 
ranging from ^ 1.5a to as large as ten times the statis- 
tical error or more. The broad range reflects the differ- 
ences from one OWLS simulation to the next. However 
the largest of these biases arc unlikely to be the product 
of any actual analysis. It seems more likely that the team 
undertaking the analysis will notice that all models pro- 
vide poor fits to the data using some "goodness-of-fit" 
criterion. Nevertheless, it remains imperative to under- 
stand the reasons for the poor fits. Our analysis suggests 
that concentration fitting may reduce systematic errors 
on dark energy equation of state parameters due to bary- 
onic effects to < 1.6(7 in the worst case and < 0.5cr in six 
of the nine simulations we have analyzed. The concomi- 
tant increase in the statistical errors is < 30%. While 
concentration fitting does alleviate biases in this case, a 
more sophisticated analysis may be necessary for data of 
the quality expected from Stage IV experiments. 

The remainder of this manuscript is organized as fol- 
lows. In the following section, we describe the lensing 
power spectra from which we aim to infer cosmological 
parameters, the details of our modeling procedure, and 
the cosmological parameters that we consider. We also 
discuss the Fisher matrix method for estimating statisti- 
cal and systematic errors in model parameters. In § IIIH 
we describe the OWLS simulations and show the differ- 
ences in lensing power spectra predicted by several of the 
simulations in the OWLS suite. We describe our simple 
mitigation model in § IIVI Our results for the statisti- 
cal and systematic errors on dark energy parameters are 
given in § |Vl where we address a DES-like experiment, 
and a future LSST- or Euclid-like experiment in turn. 
We summarize our results and present our conclusions in 

§IYI1 



II. PRELIMINARIES 
A. Weak Lensing Observables 

We consider cosmological parameter inference using 
measurements of cosmic shear from large-scale imag- 
ing surveys. We assume that each galaxy has a well- 
characterized photometric redshift estimate, so that the 
source galaxies can be binned in photometric red- 
shift bins. We infer cosmological parameters from the 
A'obs = N^iNz + l)/2 number density-weighted angular 
convergence power spectra and cross spectra among the 
galaxies in each of the redshift bins. 



three-dimensional matter power spectrum at wavenum- 
ber k and redshift z, Wi{z) are the lensing weight 
functions, and the lower-case Latin indices indicate the 
redshift bins (e.g., index i runs from 1 to iVz). The lensing 
weight functions are 

W,{z) = lnMH'oii + z)D^{z) I ^, (2) 

where dnj/dz is the redshift distribution of source galax- 
ies in the i"^ redshift bin, Hq is the present Hubble rate, 
and Da{z,z') designates the angular diameter distance 
between redshifts z and z' . 

The observed spectra P^^ , consist of terms due to signal 
(P^j) and noise, 

P^(£) = py(£) + ni5ij(7^), (3) 

where rij is the surface density of source galaxies in red- 
shift bin i, (7^) is the intrinsic source galaxy shape noise 
for each shear component, and is the Kronecker delta 
symbol. The covariance among observables is 

c[P'j{e),p^\e)] = p'^pi' + p'^pt, (4) 

assuming Gaussian statistics. Over the range of scales we 
consider, the Gaussian approximation is reasonable (e.g., 
Ref. [l^) and greatly simplifies the analysis. Moreover, 
it is a conservative assumption for our purposes because 
adopting non-Gaussian covariance generally renders sta- 
tistical errors larger and diminishes the relative impor- 
tance of the systematic errors we consider. Throughout 
this study, we adhere to a common convention by taking 

B. Survey Characteristics and Cosmological 
Parameters 

We consider cosmological constraints from two repre- 
sentative surveys. The first experiment we consider is 
based on the specifications of the Dark Energy Survey 
(DES)'^. DES is an example of a near-term, "Stage HI" 
project that will exploit cosmic shear measurements to 
derive constraints on dark energy parameters. We model 
DES by taking a fractional sky coverage of /sky = 0.12, 
corresponding to approximately 5000 deg. , and a total 
surface density of imaged galaxies of A^^ = 15 arcmin"^. 
This choice is optimistic, but it is a conservative assump- 
tion for our purposes as smaller statistical error bars set 
a more stringent requirement for the mitigation of sys- 
tematic errors. We take the DES redshift distribution 
of source galaxies from the DES Blind Cosmology Chal- 
lenge (BCC) simulation. The DES BCC simulation com- 
prises 5000 square degrees of simulated shear maps and 



In Eq.[Tl H{z) is the Hubble expansion rate, Da{z) is the 
angular diameter distance to redshift z, Ps{k,z) is the 



^ |http : //darkenergysurvey ■ org 
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is tuned to match the expected observational character- 
istics of the DES mission. We divide the source galaxies 
into Nz — 5 redshift bins such that 20% of the total num- 
ber of observed galaxies are placed in each bin, this gives 
-^obs = 15 distinct convergence spectra. Binning more 
finely in redshift does not alter our results, in accord 
with prior studies [s^l ■ Throughout the remainder of 
this paper, we will refer to results based on these survey 
specifications by the name "DES." 

In addition to DES, we will estimate the potential in- 
fiucnces of baryonic physics on dark energy constraints 
from long-term future experiments, categorized as "Stage 
IV" experiments in the report of the Dark Energy Task 
Force [23 ■ Examples of potential Stage IV experiments 
that will explore cosmological constraints from weak 
gravitational Icnsing arc the Large Synoptic Survey Tele- 
scope (LSST, Ref. ipil')^ or the European Space Agency's 
Euclid^ project [321 ■ We characterize these experiments 
by a fractional sky coverage of /sky = 0.5 and a number 
density of source galaxies of N-^ = 30 arcmin^^ ^. Again, 
these choices are optimistic, but they maximize the rel- 
ative importance of the systematics we aim to militate 
against, so they are conservative choices for our purposes. 
We model the redshift distribution of source galaxies in 
these long-term surveys as dn/dz oc exp(— (z/zq)^'^), 
with zq ~ 0.34 to give a median redshift Zmcdian ~ 1- 
This choice is based on the approximate, observed dis- 
tribution of high- redshift galaxies [s^l- As with DES, we 
place the source galaxies into = 5 redshift bins so 
that the 20% of the galaxies fall into each bin. We refer 
to results with these specifications as "Stage IV" results. 

We consider cosmologies defined by seven parameters, 
three of which describe the dark energy. The parameters 
that describe the dark energy are the contemporary dark 
energy density in units of the critical density, flDE, and 
the two parameters wq and that specify a dark energy 
equation of state that varies linearly with scale factor, 
u'(a) = wo+Wi^{l—a). The parameters of our fiducial cos- 
mology are fixed to match the cosmological parameters 
assumed in the OWLS simulation program [26| . The pa- 
rameters specified in the OWLS program are the matter 
density, ujm = 0.1268, the baryon density, ujb = 0.0223, 
the scalar spectral index, rig = 0.951, the amplitude of 
curvature fluctuations on a scale of /c = 0.05 Mpc~^, 
= 1.9 X 10~^ (we actually vary InA^ about this 
value), Ode = 0.762, wq = —1, and = 0. These pa- 
rameter values imply that the root-mean-square matter 
density fluctuation on a scale of 8 Mpc is trg = 0.74. 
We include prior constraints on these parameters that re- 
flect expected limits from the Planck cosmic microwave 
background anisotropy measurements in all of our calcu- 



* http : //www. Isst ■ org I 

^ http : //sci . esa. int/euclid| 

® The sky coverage of Euclid will likely be closer to /^j^y « 1/3 



lations. The Planck prior matrix that we use was com- 
puted in Ref. [13] ■ In addition to these seven cosmo- 
logical parameters, we introduce three other parameters, 
described in § llVi that account for baryonic effects. 



C. Methodology 

In principle, we propose to assess the effectiveness of 
the mitigation approach proposed in Ref. [2l[ using the 
following steps. 

1. Take the lensing spectra predicted by one of the 
OWLS simulations as mock data. 

2. Fit the mock data to a model by varying 7 cos- 
mological parameters and determine the statistical 
errors on the cosmological parameters. This model 
does not include the effects of baryons. 

3. Determine the raw bias as the difference between 
the resulting best-fit dark energy parameters and 
the "true" input parameters used to generate the 
OWLS simulations. 

4. Fit the mock data again to a model with those same 
7 cosmological parameters as well as 3 additional 
parameters that account for baryonic effects. 

5. Determine the residual bias as the difference be- 
tween this second fit and the "true" values of the 
dark energy parameters used in the OWLS simula- 
tions. 

6. Compare the size of the error bars in both cases to 
see the amount by which the errors are inflated as 
a result of the new degrees of freedom. 

7. Repeat for each of the OWLS simulations to arrive 
at nine distinct assessments. 

In practice, going through this entire process for all the 
cases of interest would be extremely time consuming, be- 
cause fitting for the cosmological and concentration pa- 
rameters in the multi-dimensional parameter space that 
we explore is a computationally-expensive task. Instead, 
we proceed using an approximation, based on both di- 
rect fitting for model parameters and Fisher matrix (de- 
scribed below) estimates for the statistical and system- 
atic errors in model parameters. However, it is important 
to stress that we are assessing the residual bias that will 
ensue if analysts follow the mitigation strategy proposed 
in Ref. [2l| on upcoming data sets. 

In order to limit computational effort, we use the 
Fisher information matrix to assess the constraining 
power of these A'obs observable spectra. Wc assume that 
the spectra are independent, Gaussian random variables 
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at each multipole, so that the Fisher matrix may be writ- 
ten as 



AB 



^max 

fc^min i=l j=i k=l l=k 



dpA 



on 

dpB 



(5) 



The matrix represents the prior constraints on the 
cosmological parameters, pa are the parameters of the 
model, and C~^[P^ , P^^] is the inverse of the covariance 
matrix between observables. The upper-case Latin in- 
dices signify model parameters. The parameter /sky spec- 
ifies the fraction of sky observed in the survey, and the 
sum runs over multipoles from to ^max- We take 

— 1/2 

imin = 2/sj.y ; however all of the constraints we con- 
sider are dominated by multipoles significantly higher 
than fmin so that this choice is inconsequential. We 
take ^max = 5000 throughout most of our study so as 
to remain in a regime in which a number of simplif ying 
assumptions are approximately valid (e.g., Refs. |29l. l35l- 
Isij). but we explore other choices of maximum multipole. 
Including such high multipoles in our analysis may well 
be overly optimistic. However, using higher multipoles 
(smaller scales) in the cosmological analysis results in 
greater constraining power, so it is interesting to deter- 
mine the utility of our mitigation scheme out to relatively 
high multipoles. The Fisher matrix approximates the co- 
variance among model parameters at the maximum of the 
likelihood, so that the error in the estimate of the A**^ pa- 
rameter can be approximated as cr(pA) — \J -^aA' ^-ftsr 
marginalizing over the remaining parameters. 

The Fisher matrix formalism provides a straightfor- 
ward estimate of parameter biases due to undiagnosed, 
systematic offsets in observables. Let AP^j represent the 
difference between the true observable and the observable 
perturbed due to some systematic effect. A Taylor ex- 
pansion about the maximum likelihood gives an estimate 
of the systematic error contribution to model parameter 
PA due to the systematic offsets in observables [40l |: 



b{pA) - E^abE(2^+1)/^ 



sky 



B 



i,j,k,l 



dpB ' 



(6) 



The sums over observable (lower-case Latin) indices in 
Eq. ([6]) have the same form as those in Eq. ([5]), though 
we have written them as a single sum for brevity. 

The practical strategy that we implement in an effort 
to limit computational expense is a modification to the 
ideal strategy that we would, in principle, pursue as de- 
scribe above. We repeal and replace Step 2 through Step 
5 with the following steps. 



2'. Determine the raw bias using the Fisher matrix re- 
lation [Eq. dH])] with the systematic offsets in power 
spectra given by the difference between the OWLS 
baryonic simulation and the fiducial OWLS simu- 
lation that treats dark matter only. 

3'. Fit the 3 parameters of the concentration-mass re- 
lation of halos to the OWLS convergence power 
spectra within a fixed cosmological model. This re- 
sults in a best-fit concentration-mass relation that 
best describes the OWLS simulations within the 
true underlying OWLS cosmological model. 

4'. Compute the residual differences between the 
predicted convergence spectra from the OWLS 
simulations and those predicted by the best-fit 
concentration-mass model with cosmological pa- 
rameters fixed to the OWLS cosmological model. 
These residual differences, AP^j, will give rise to 
systematic errors in inferred cosmological parame- 
ters. 

5'. Use the Fisher matrix formalism [particularly 
Eq. ([6])] to estimate the residual bias from AP^^ 
after allowing for parametric freedom in the 
concentration-mass relation. This level of bias will 
not be removed by fitting for halo structure and 
will remain in the error budget on cosmological pa- 
rameters. 

This strategy enables us to estimate the efficacy of the 
baryonic physics mitigation proposal. Additionally, it 
requires explicitly fitting for only the concentration-mass 
relation (a 3-dimensional subspace of the full parame- 
ter space), so that the computational effort required is 
significantly less than would be required to fit for halo 
structure and cosmology simultaneously. The merit of 
this approach is that the fit allows us to assess the fi- 
delity with which we can model lensing power spectra in- 
cluding baryonic effects and it brings our fiducial model 
closer to the true, underlying model prior to applying the 
Fisher formalism in step 5'. This modified procedure in- 
troduces a possible source of confusion. One might think 
that our final statistical and systematic error estimates 
do not include the covariance between cosmological and 
concentration parameters. We emphasize that this is not 
the case. The Fisher matrix formalism for computing 
biases in model parameters accounts for the covariance 
among all the model parameters. The limitation is the 
standard caveat that the Fisher matrix is a first-order 
approximation about the maximum likelihood. 



III. SPECTRA FROM SIMULATIONS 

We estimate the infiuencc of baryonic processes on 
weak lensing power spectra using the Overwhelmingly 
Large Simulations (OWLS) [13, HI, [13 . The OWLS suite 
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includes a set of ten simulations set within the same cos- 
mological model that evolve from the same initial con- 
ditions to redshift z = 0. One of these simulations, the 
"DMONLY" simulation, treats the gravitational evolu- 
tion of dissipationless dark matter only. This type of dis- 
sipationless dark matter simulation is similar to the vast 
majority of simulations that have been used to model the 
observations of contemporary and future surveys. The 
remaining nine simulations include the baryonic compo- 
nent of the Universe along with various effective models 
for baryonic gas cooling, star formation, and a number 
of feedback processes. It is not feasible to simulate these 
processes directly, so baryonic simulations rely on a vari- 
ety of effective models for these processes. Effective mod- 
els for baryonic processes remain quite uncertain and it 
is not possible to produce a definitive prediction for the 
influences of baryonic processes on observables, such as 
convergence power spectra. The utility of a simulation 
suite such as OWLS, is that it provides a range of dis- 
tinct, but plausible predictions for observables so that 
systematic errors induced by our ignorance of baryonic 
physics can be estimated. An important advantage of the 
test that we present here is that we are applying a mit- 
igation strategy developed on the simulations of Rudd 
et al. [3 to an independent set of simulations that were 
performed using markedly different simulation strategies. 
The details of the OWLS simulations have been given in 
Refs. [l^, [4l| - |45j . while the OWLS power spectra were 
the subject of Ref. [l^l, to which we refer the reader 
as these details are not of immediate importance in the 
present paper. Follow up studies by McCarthy et al. [i^ 
and McCarthy et al. [43] suggest that the properties of 
galaxies and hot gas in galaxy groups are modeled most 
reliably in the "AGN" simulation, which includes strong 
feedback from active galactic nuclei. 

We regard each of the nine baryonic simulations sepa- 
rately as a potential realization of the effects of baryons 
on convergence power spectra. Accordingly, we treat the 
convergence spectra from each of the OWLS baryonic 
simulations as "observed" spectra that must be modeled 
faithfully in order to extract reliable constraints. The aim 
is to test whether a specific strategy to mitigate the in- 
fluence of baryonic processes on dark energy constraints 
can be applied to a variety of distinct predictions success- 
fully. Success in this context means that the mitigation 
procedure renders the biases in dark energy parameters 
significantly smaller than the expected statistical errors. 
If a single mitigation strategy were to achieve the req- 
uisite reduction in dark energy parameter biases for all 
plausible simulations, it would be a strong indication that 
the mitigation strategy may be applied to observational 
data to limit systematic errors on dark energy parameters 
associated with the influences of baryons. 

van Daalen et al. [l^l used the OWLS simulations 
to study the effects of baryonic physics on the matter 
power spectrum. We use the 3D matter power spec- 
tra, Ps{k,z), provided in Ref. [13] for the OWLS sim- 
ulations to estimate convergence power spectra using 



Eq. (IT|). In practice, the tabulated matter power spectra 
from Ref. [23| cannot be used directly to predict conver- 
gence power spectra. Due to computational limitations, 
the OWLS simulation volumes are relatively small (cu- 
bic boxes L = 100 h^^Mpc on a side), and are subject to 
significant sample variance and finite volume effects on 
large scales. In order to overcome these drawbacks, we 
utilize the OWLS Ps{k,z) tables directly for wavenum- 
bers k > 0.314 hMpc~^. The OWLS spectra are reliable 
for k < 10 /iMpc~^ [131, which is sufficient for our pur- 
poses [l3|- For wavenumbers k < 0.314, we use the halo 
model as implemented in Ref. [ll| to estimate the mat- 
ter power spectrum. We multiply the halo model power 
spectra by a correction factor that ensures that the two 
spectra agree a,t k = 0.314 hMpc~^ . In the OWLS simu- 
lations, baryonic effects induce changes in power spectra 
of < 1% on scales k < 0.314 /iMpc~^. An important 
caveat to our approach is that assume that the effects of 
baryons on scales k < 0.3 hMpc~^ are insignificant. 

Figure [5] shows the fractional differences between the 
convergence power spectra predicted by the baryonic sim- 
ulations compared to the DMONLY simulation. For sim- 
plicity, we show only the power spectrum in our third 
DES redshift bin, P^a (0.5025 <z< 0.6725 for our DES 
model). The predictions for the other fourteen observ- 
ables show similar features. Fig. [2] shows these residuals 
in several distinct ways. The shaded band is the envelope 
of the power spectrum residual constructed from all nine 
of the baryonic simulations. The residual power spectra 
from three specific simulations, namely AGN (bounding 
the shaded region above), NOSN (bounding the shaded 
region below for £ > 2000), and WDENS (intermediate), 
are shown as solid lines. 

As is evident in the AGN, NOSN, and WDENS 
models, the deviations in the spectra predicted by 
the baryonic simulations differ from the DMONLY 
simulation in a way that is correlated from multipole 
to multipole. Accordingly, the shaded bands in Fig. [2] 
represent the envelope of deviations, while an individual 
spectrum will not range over the shaded To 
represent the typical shapes of the baryonic simulation 
spectra, we have performed the following exercise. 
We have computed the covariance among the distinct 
spectra of the baryonic simulations, C[Pj^^{£), Pj^'^{£')] = 

Ell (^f ''w - {p'.'m {p'.''Kn - (p'jin)), 

where i is an index designating the simulation, iV = 9 
is the number of baryonic simulations, and {Pi^^{£)) 
is the average of the spectra from all of the simula- 
tions. We then diagonalized the covariance matrix, 
C[P^'^{£), Pj^^ {£')]. The eigenvectors of the covariance 
matrix represent the principal modes of variation of the 
power spectra and the eigenvalues represent the variance 
accounted for by the corresponding eigenvectors. The 
dashed lines in Fig. [5] are the two eigenvectors corre- 
sponding to the first- and second-largest eigenvalues. 
These principal modes account for over 90% of the 
variance among the spectra and illustrate the correlated 
manner in which the baryonic simulation spectra may 
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FIG. 2: Convergence power spectrum residuals between 
the DMONLY simulation and baryonic simulations (denoted 
"BAR" in the vertical axis label). For simplicity, this fig- 
ure shows only the residual for the auto power spectra in the 
third tomographic bin (0.5025 < z < 0.6725), P^^, assum- 
ing a DES analysis. Fractional residuals of the other fourteen 
observables display similar features. The shaded band cov- 
ers the region spanned by the residuals of all of the OWLS 
baryonic simulations. The three thick, solid lines show three 
specific simulations that contribute to this band, namely, the 
AGN (bounding the shaded region above), NOSN (bounding 
the shaded region below for £ > 2000, and WDENS (interme- 
diate) simulations. The residuals at different multipoles are 
highly correlated. The dashed lines show the principal modes 
of the residuals that have the highest (upper) and second- 
highest (lower) variance. These modes account for over 90% 
of the variance among the spectra and demonstrate the cor- 
related manner in which baryons alter lensing power spectra. 



FIG. 3: Convergence power spectrum residuals between our 
halo model fits to the spectra predicted by the baryonic sim- 
ulations and the baryonic simulation predictions themselves. 
The shaded band shows the envelope containing the fit residu- 
als for all of the OWLS baryonic simulations, analogous to the 
shaded band in Fig. O As in Fig. [21 the thick, solid lines show 
examples of the fit residuals from specific cases, namely the 
AGN (lowest at high multipole), NOSNJ^OZCOOL (highest 
at high muhipole, lowest at ^ < 1000), and WDENS (interme- 
diate) simulations. It is evident that the envelope containing 
all residuals is comprised by different simulations over differ- 
ent multipole ranges. The dashed lines show the two principle 
modes of the residuals with the highest variance. These rep- 
resent correlated failure modes of the fitting procedure. This 
figure is to be compared to Fig. [2] but notice that the ranges 
on the ordinal axes are considerably different in the two fig- 
ures. The halo model fits typically reduce power spectrum 
residuals by more than an order of magnitude. 



differ from the DMONLY predictions. 



IV. FITTING FOR BARYONIC EFFECTS WITH 
HALO CONCENTRATIONS 

Motivated by prior studies indicating that the largest 
effect of baryons on convergence spectra on relevant 
scales is a modification of halo structure [l^, [2l| , we pur- 
sue a mitigation strategy in which baryonic effects are 
entirely encapsulated into changes in the internal mass 
distributions within dark matter halos. It is certainly 
not true that the only effect of baryons is to alter halo 
structures. For example, the distribution of halo masses 
changes slightly (e.g., Ref. [H, El]), and baryonic effects 
extend beyond halo virial radii (e.g., Ref. Our 



goal is to determine the practical utility of such a model 
in analyses of forthcoming data. 

We assume that the average mass distributions within 
dark matter halos can be described by the density profile 
of Ref. m (NFW hereafter). 



' (Cr/i?200m)(l+cr/i?200m)'' 

The parameter c is the halo concentration and the den- 
sity profile is normalized by our definition of a halo as 
a spherical object within which the mean density is 200 
times the mean density of the universe, pvu- = 200/9m- 
Therefore, halo mass and radius are related by m = 
47rpvirJ^200m/'^i that the profile can be normalized by 
m = An /g^^""" p{r) dr for a given mass and concen- 
tration. Halo concentrations predicted by dissipation- 
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TABLE I: Concentration parameters that best-fit the OWLS 
simulation spectra. The best-fit values are given for the nor- 
malization of the concentration relation co, in units of the 
normalization for the DMONLY simulation, the power-law 
index specifying the mass dependence a, and the power-law 
index of the redshift dependence, j3. The values in this table 
correspond to the specific case of DES with £max = 5000. The 
other cases were explore yield very similar best-fit concentra- 
tions. 



of iVz convergence power spectra. We fit the spectra by 
minimizing 



Simulation 


/ DMDNLY 
CO /Co 


a P 


REF 


1.13 


0.20 0.83 


WML4 


1.14 


0.18 0.91 


NOSN 


1.20 


0.27 0.50 


NOZCOOL 


1.16 


0.16 0.91 


WDENS 


0.96 


0.20 1.39 


WML1V848 


1.09 


0.22 1.33 


DBLIMFV1618 


0.89 


0.24 1.23 


NOSNJIOZCOOL 


1.41 


0.40 1.44 


AGN 


0.34 


0.79 2.06 



less simulations of dark matter only have been stud- 
ied extensively (recent examples are Refs. [13, [1^ IHol - 
[5^). The relationship between halo concentration, halo 
mass, and redshift in the OWLS DMONLY simulation 
can be adequately characterized by a power-law distri- 
bution [23,[5i,[5|^ 



M 

ciM,z) = co[ — 



(1 + ^) 



-13 



(8) 



with Co = 7.5, a = 0.08, and /3 = 1, in broad agreement 
with prior studies. The parameter Mp is a pivot mass, 
which we take to be Mp = 8 X 10^^ /i-^Mq. We choose 
the pivot mass to be close to the halo mass that is most 
well constrained by lensing spectra [21 1. 

We describe modifications to halo structure through a 
modified concentration relation, following Refs. [l8ll2l|. 
In particular, we allow the parameters cg, a, and /? in 
Eq. ([8|) to vary in order to describe convergence power 
spectra within the baryonic OWLS simulations. We de- 
termine the values that best capture the simulation re- 
sults as follows. For each simulation, we produce a set 



The fit values come from Ref. [ssll after applying a correction to 
change the pivot mass to Mp at 2 = 0. Ref. [53| explored a set of 
simulations that differ slightly from the OWLS simulations that 
we examine here. Ref. [54j and Ref. [27j confirm that the OWLS 
DMONLY simulation has a very similar concentration-mass re- 
lation, but do not provide detailed fits to the concentration-mass 
relation as a function of redshift. For the purposes of the present 
paper, this relation serves only to establish a baseline with re- 
spect to which we model the remaining baryonic simulations. 
Plausible alterations to the baseline model, given the results 
quoted in Refs. [53L l54|. cause only minor quantitative changes 
to our subsequent results. 



= £ (2£+l)/.ky 

^[P^(£),Pf(£)]<5Pf(^), (9) 



i,j,k,l 



where SP^ {£) is the difference between the model and the 
simulation prediction, for the concentration parameters 
Co, a, and /? at fixed cosmology. This results in best- 
fit values for the concentration parameters, and residual 
differences between the best-fit modified concentration 
models and the predicted spectra from the baryonic sim- 
ulations. In this manner, we assess the ability of the 
modified concentration model to describe the baryonic 
simulations. 

The implementation of our model for the effect of 
modified concentrations on convergence power spectra 
is based upon the halo model for Ps{k,z). The de- 
tails of the halo model implementation are described in 
ReL 



21 



However, using the halo model to fit con- 
vergence power spectra from the OWLS simulations in- 
troduces a non-trivial complication. On large scales, 
the halo model predictions for Ps{k,z) are systemati- 
cally offset from the OWLS simulation predictions be- 
cause of the finite volumes of the simulations. To over- 
come this, we proceed as follows. For each trial value of 
the halo concentration parameters, we compute a mat- 
ter power spectrum offset APs^nuik, z) between the fidu- 
cial halo model with the standard values of co, a, and 
(3 and the halo model prediction with our trial values 
of the concentration parameters. We then compute our 
trial matter power spectrum, which we use to predict 
convergence and compute x^, by adding the offset de- 
fined by the halo model to the DMONLY prediction, 
P5,trial(fc, z) = P5,DMONLY(fc, z) + APs^Huik, z). In Other 
words, we utilize the halo model to estimate a correction 
to be applied to the DMONLY matter power spectra. 
In this manner, the spectra are not offset systematically 
with respect to each other on large scales. This strategy 
mimics what would likely have to be done in any anal- 
ysis of real data; the predictions of dissipationless dark 
matter simulations would be established, but a correc- 
tion would need to be applied to account for baryonic 
effects. This is the same strategy suggested in Ref. [l^ 
and adopted by Ref. 

Figure [3] summarizes the results we obtained by fit- 
ting concentrations to describe the baryonic OWLS sim- 
ulations. The figure shows the difference between the 
power spectra in the best-fit models with modified con- 
centrations and the "true" baryonic simulations. This 
is analogous to Fig. [51 As in Fig. [51 Fig. [31 shows only 
residuals of the auto-correlation for sources in the third 
DES photometric redshift bin, in the interest of clarity. 
The most obvious feature of Fig. [31 is that the fit resid- 
uals are nearly an order of magnitude smaller than the 
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differences between the baryonic and DMONLY simula- 
tions. The correlated structure of the fit residuals is also 
shown in Fig. |3]by the two principal modes that account 
for the majority of the variance in the spectra among the 
baryonic OWLS simulations. The two modes depicted in 
Fig. [3] account for nearly 97% of the variance among the 
spectra in the third rcdshift bin. 

The residuals in Fig. [3] represent the remaining sys- 
tematic errors in the predictions of convergence power 
spectra after accounting for baryonic effects using a sim- 
ple, phenomenological model. These systematic offsets 
will propagate into biases in the estimators of the cosmo- 
logical parameters. We explore this in detail in the next 
section. 

The best-fit concentrations are given in Table H] The 
concentrations of the dark matter halos in the OWLS 
simulations have been presented in Duffy ct al. (H^j . The 
concentrations we quote in Table |T] exhibit two important 
deviations from the results quoted in Ref. [H^]. First, 
our fit to the AGN simulation results in significantly 
lower concentrations. This could be because mass that 
would be within the DMONLY virial radii of these halos 
is moved outside of the virial radii in the AGN simula- 
tion and the concentrations are driven to low values to 
account for this effect. The OWLS collaboration has ar- 
gued that the AGN simulation is the most realistic among 
their suite, so a genuine discrepancy in this case could 
have important implications for our method. Second, 
our fits generally yield a stronger mass dependence than 
the OWLS fits. 

A comparison of between our concentration results is 
not trivial for several reasons. Foremost among these, 
Ref. [13] quote the effective concentrations of dark mat- 
ter while Icnsing is sensitive to all matter. The concen- 
trations that we quote are the effective concentrations of 
all matter (baryonic as well as dark). Ref. [13] does pro- 
vide non-parametric measures of total matter concentra- 
tion; however these measures exhibit behavior different 
from the NFW fits and it is unclear how they correspond 
to the result of our fitting procedure. Additionally, the 
lensing signal on the scales that we explore is sensitive 
to mass redistribution on scales near halo virial radii and 
insensitive to halo profiles on significantly smaller scales. 
Therefore, it is not clear that our exercise should yield 
the same concentrations as those derived from fitting pro- 
files directly because the two procedures are not equally 
sensitive to halo structure at all scales. Furthermore, 
we have assumed power law relationships between halo 
concentration, halo mass, and rcdshift. This can be jus- 
tified by the fact that Icnsing is sensitive to a relatively 
narrow range halo masses and redshifts [2l| , so that the 
power-law indices that we recover may not correspond to 
those derived from a fit to simulation results over a wide 
range of masses and redshifts. As a results of these com- 
plicating factors, we reserve a more detailed comparison 
between our concentration parameters and those of the 
simulations for future work. 



V. COSMOLOGICAL CONSTRAINTS AND 
RESIDUAL BIAS 

In this section, we project the effects of baryons in 
simulations onto cosmological parameters. We begin by 
discussing our results in the context of the DES and con- 
clude the section with a brief discussion of possible cos- 
mological biases for stage IV experiments such as LSST 
or Euchd. 

It is useful to consider the baseline constraints on the 
dark energy equation of state parameters. We consider 
two cases that will prove useful in the following. The 
first case corresponds to standard constraint projections 
on the dark energy equation of state parameters assum- 
ing that the nonlinear growth of structure is known per- 
fectly. In the context of our analysis, this means that the 
concentrations of halos are known perfectly and we re- 
fer to these constraints as "Fixed-C" constraints accord- 
ingly. As we are exploring a mitigation strategy in which 
we fit for concentrations concurrently with cosmological 
parameters, it is necessary to assess the degradation in 
dark energy parameter constraints due to this additional 
freedom. We refer to constraints derived from an analy- 
sis in which concentrations are fit alongside cosmological 
parameters as "Fit-C" . 



A. Results for the Dark Energy Survey 

Our baseline constraints for DES, as well as the degra- 
dation in constraints when concentrations are fit along- 
side cosmological parameters, are shown as a function of 
the maximum multipole used in the analysis in Figure SI 
Fig. |4] contains four panels. The two panels on the left 
show results for wq, while the two panels on the right 
show results for Wn. The top panels show the marginal- 
ized constraints on the equation of state parameters as 
a function of the maximum multipole used in the analy- 
sis. The solid lines show the marginalized constraints in 
the standard Fixed-C case, while the dashed lines show 
constraints in the Fit-C case in which concentrations are 
permitted to vary. The lower panels show the ratio of 
the Fit-C constraints to the Fixed-C constraints at each 
multipole, giving the factor by which introduction of the 
additional nuisance parameters describing concentrations 
degrades the constraints. 

It is clear that constraints are degraded if the 
concentration-mass relation of dark matter halos must 
be allowed to vary. This degradation is mild {< 20%) if 
the maximum multipole used in the cosmological analy- 
sis is ^max ^ 3000, and increases to > 40% once scales 
to £max ~ 5000 are included. The level of degradation 
depicted in Fig. 2] is slightly larger than that estimated 
in Ref. [2l| for a DES-like experiment. We find that 
this discrepancy is almost entirely due to the fact that 
the source redshift distribution used in the present study 
differs significantly from that assumed in Ref. [2l|. In 
particular, the redshift distribution that we assume con- 
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FIG. 4: Constraints as a function of maximum multipole used to infer cosmological parameters. The left panels show results 
for wo, while the right panels show results for w^. In the top panels, the solid lines show standard constraints assuming that 
the power spectrum is known perfectly. The dashed lines show the constraints that may be realized if the halo concentration- 
mass relation must be fit simultaneously with cosmological parameters. The bottom panels show the ratio of the constraints 
realized when concentrations must be fit, to the constraints realized in the standard scenario, cr^"""*~*-'. Constraints are 

significantly degraded beyond ^max ^ 3000 when the concentration-mass relation is permitted to vary. Constraints on wo and 
TOa show the same qualitative features. 



centrates source galaxies at significantly lower redshift, 
resulting in relatively lower lensing power compared to 
noise and reducing the lever arm to high redshift sources. 

To determine the impact of this mitigation scheme 
on the dark energy program of a particular experiment, 
such as DES, the degradation in cosmological param- 
eters caused by fitting for the concentrations of halos 
must be compared to the biases in these parameters that 
may be realized if baryonic effects are neglected. We 
compute these biases by using the residuals between the 
DMONLY and baryonic simulations as the systematic 
offsets in Eq. (jH]). An example of these offsets is depicted 
in Fig. l 

The maximal biases that arc induced by neglecting 



baryonic effects in our analysis of the OWLS simulations 
are shown as the outer (blue) bands in Fig. [S] (for wq) 
and Fig. |6] (for Wa). In particular, the outer {blue) bands 
delineate the extremal biases (maximum and minimum 
as the biases may be positive or negative) induced by an- 
alyzing any of the OWLS simulations without account- 
ing for baryonic effects. The outer bands in Fig. [5] and 
Fig. ini represent the envelope of the bias from all simu- 
lations, while for any individual simulation, the bias is a 
smooth function of maximum multipole. The features in 
the bands in Fig. [5] and Fig. |6] arise when the particular 
simulation that gives rise to the extremal bias changes 
from one multipole to the next. 

The biases induced by neglecting baryonic effects in 
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FIG. 5: Biases induced in the dark energy equation of state 
parameter wq as a function of the maximum multipole used 
to infer cosmological parameters. The bias is sliown in units 
of the statistical error on wq in order to make the relative 
importance of the systematic error induced by baryons ex- 
phcit. The outer shaded (blue) band covers the range of biases 
spanned by computing the biases induced by analyzing all of 
the OWLS baryonic simulations without any model for bary- 
onic effects. The inner shaded (orange) band shows the range 
of biases induced by after taking the best-fit concentration 
model to describe baryonic effects in the OWLS simulations. 



FIG. 6: Biases induced in the dark energy equation of state 
parameter uia as a function of the maximum multipole used 
to infer cosmological parameters. The bias is shown in units 
of the statistical error on Wa in order to make the relative 
importance of the systematic error induced by baryons ex- 
plicit. The outer shaded (blue) band covers the range of biases 
spanned by computing the biases induced by analyzing all of 
the OWLS baryonic simulations without any model for bary- 
onic effects. The inner shaded (orange) band shows the range 
of biases induced by after taking the best-fit concentration 
model to describe baryonic effects in the OWLS simulations. 



Fig. [5] and Fig. IH] arc as large as Sc for wq and Wa, 
if the analysis include multipoles out to ^max ~ 3000. 
Including multipoles out to i^max ~ 5000 drives the max- 
imal potential bias to 6a. These bias levels have a 
clear significance for the effort to understand dark en- 
ergy. However, we remind the reader that we have used 
the Fisher matrix approximation to estimate the biases 
on cosmological parameters. A necessary caveat to our 
results is that Eq. ^ is the lowest-order approximation 
to the bias in the limit of small parameter biases and 
may not provide an accurate bias estimate for large bi- 
ases. Nevertheless, the statement that the biases arc sig- 
nificant (> la) in this case is robust. 

The inner (orange) bands in Fig. [5] and Fig. [B] de- 
lineate the extremal range of dark energy equation of 
state parameter biases realized after fitting for the halo 
concentration-mass relation in the baryonic simulations 
and using the fit to correct the power spectra as described 
in § HyI Specifically, we compute these biases by utihz- 
ing the residuals between the corrected DMONLY spec- 
tra and the baryonic simulations, an example of which is 
shown in Fig. |31 in Eq. (jH]). To make Fig. [5] and Fig. [5] 
show fair comparisons of the biases, the biases for the fit- 



corrected cases are shown in units of the statistical error 
in the case of fixed concentrations ("Fixed-C"). Showing 
these biases in units of the statistical error in the case 
of varying concentrations ("Fit-C") would reduce their 
magnitudes in Fig. [5] and Fig. [51 

Fitting for concentrations clearly leads to dramatic re- 
ductions in parameter biases. Indeed, the biases are typi- 
cally less than ~ 10% of the statistical error at low multi- 
poles and never exceed ~ 50% 60%) of the statistical 
error in wq (wa) for ^max 5000. This suggests that the 
mitigation strategy of fitting for a halo concentration re- 
lation alongside cosmological parameters will result in a 
dark energy error budget that is preferable to neglecting 
baryonic effects. At .^max ^ 3000, fitting for concentra- 
tion increases the statistical error in wq, for example, by 
~ 20% compared to the ideal case (Fig. |4]) , and reduces 
the systematic error to ~ 40% of the statistical error 
(maximum) . Taking a simple and conservative approach 
of adding these two contributions, the resulting error on 
Wo increases to ~ 160% of the constraint in the ideal case. 
This is to be compared to a potential systematic error in 
the case where no mitigation for baryonic processes is 
undertaken of as much as ^ 350% of the statistical error. 
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B. Stage IV Experiments 

Having discussed the utility of the proposed mitigation 
scheme of Zentner et al. [llj for DES, we briefly describe 
analogous results for forthcoming, Stage IV dark energy 
experiments such as LSST and Euclid. Figure [7] shows 
cosmological constraints for a Stage IV dark energy ex- 
periment, including the degradation in those constraints 
incurred by fitting for the concentrations of halos concur- 
rently with the cosmological parameters. Fig. [7] exhibits 
two notable features compared to the analogous results 
for DES (Fig.S]). First, the constraints from Stage IV ex- 
periments are significantly more restrictive, though this 
is an expected result (e.g., Ref. [ll])- Second, Stage IV 
experiments suffer from slightly greater degradation in 
dark energy parameter constraints when fitting for halo 
structure along with cosmological parameters. In partic- 
ular, the degradation in dark energy parameters reaches 
~ 30% at a maximum multipole of ^max ~ 3000 and 
w 55% at ^,„ax ~ 5000. 

The biases resulting from the analysis of STAGE IV 
experiments are shown in Figure [8] (wq) and Fig. |9](wa). 
Notice that we present the results in a slightly differ- 
ent way for STAGE IV experiments. In particular, four 
of the OWLS simulations ("AGN", "DBLIMFV1618" , 
"WDENS", and "WML1V848") give significantly larger 
biases than any of the other five simulations. Therefore, 
we show the biases for these individually in the main 
panels of Fig. [8] and Fig. [9l We show results for the re- 
maining five simulations in the inset panels. In all cases, 
the value of fitting concentrations to mitigate for bary- 
onic effects in cosmological parameter analyses is appar- 
ent. However, notice that the residual biases in the worst 
cases can remain significant compared to the ideal sta- 
tistical error even after fitting for the concentration-mass 
relation. This indicates that a more accurate mitigation 
scheme will be needed in order to reduce the theoretical 
systematic associated with baryonic physics to the level 
of the statistical errors expected of STAGE IV dark en- 
ergy experiments. 

VI. CONCLUSIONS 

We have explored the viability of a strategy to miti- 
gate the influence of baryonic effects on dark energy con- 
straints from cosmological weak Icnsing. The strategy 
entails fitting lensing data for both cosmological param- 
eters and the concentration-mass relation of all matter in 
halos simultaneously. We assessed this scheme by using 
it to analyze power spectra predicted by a suite of cos- 
mological simulations as though they were genuine data. 
Specifically, we computed the resultant systematic and 
statistical errors on dark energy parameters that would 
be incurred by fitting such data according to this strat- 
egy- 

We find that introducing additional parameter freedom 
to describe the concentration-mass relations of halos re- 



duces the systematic errors on dark energy parameters 
to marginally-acceptable levels in all cases. For a DES- 
like analysis exploiting all multipoles £ < 3000 (5000), 
fitting for concentrations increases statistical errors by 
< 25% (45%) (Fig.g]), but reduces the potential system- 
atic error by a factor of as much as ~ 7, to < 0.3cr (O.Scr), 
in the worst case scenario (Fig. [5] and Fig. [51). The reduc- 
tion in systematic error outweighs the increase in statis- 
tical error suggesting that this mitigation scheme may be 
a viable option for analyzing data of the quality expected 
from DES. 

For Stage IV experiments, such as the surveys to be un- 
dertaken by the LSST or Euclid, the conclusion is some- 
what less straightforward. What is clear is that some mit- 
igation strategy for baryonic effects is necessary. System- 
atic errors on wq and Wa incurred by fitting the OWLS 
simulations with no model for baryonic processes can be 
as large as several tens of the statistical error if all scales 
£ < 5000 are included in the analysis (Fig. [S] and Fig. [51). 
Of course, as we mentioned in the introductory section, 
such large biases are not likely to be realized as the result 
of any analysis. Rather, it is likely that the analysis team 
will not find acceptable fits to the observables according 
to a specific fit criterion. Nevertheless, it is apparent that 
a model for possible baryonic effects will be necessary in 
order to extract cosmological parameters reliably. 

For Stage IV experiments, the strategy of fitting the 
concentrations of halos in order to militate against large 
biases in the inferred cosmological parameters, particu- 
larly the dark energy equation of state parameters, is 
relatively less effective. One complicating factor is that 
the different OWLS simulations lead to more disparate 
conclusions in this case. In the worst case, that of the 
OWLS "AGN" simulation, the residual biases after fit- 
ting for concentrations are ^ 1.6(7, assuming all scales to 
^max = 5000 are included. It is necessary to restrict con- 
sideration to multipoles inmx ^ 1100 in order to reduce 
this bias to ~ la. However, for six of the nine OWLS sim- 
ulations that we have analyzed, the residual bias includ- 
ing all scales to ^max = 5000 is < 0.5cr. The concomitant 
cost of the additional parameters for the statistical errors 
is < 55%. Fitting for an effective halo concentration- 
mass relation does reduce biases in the dark energy equa- 
tion of state parameters; however, in the most extreme 
cases that we have analyzed, these biases remain signif- 
icant compared to the statistical errors expected from 
Stage IV experiments. 

As pointed out by Zentner ct al. [2l|, fitting for an 
effective concentration-mass relation also yields informa- 
tion that may help to constrain galaxy formation models. 
In this case, the procedure gives constraints on the pa- 
rameters of the concentration-mass relation at no addi- 
tional cost. In the case of DES (Stage IV), the best- 
constrained halos have masses M ^ 8 x 10^'^/i~^Mo 
(Af ~ 6 X lO^^/i-^Mo) at redshift z ~ 0.23 (z - 0.31) 
and constraints on the average concentrations of such ha- 
los are ac/c ~ 0.06 [ocjc ~ 0.03). Such constraints may 
prove useful in understanding the formation histories of 
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FIG. 7: Constraints as a function of maximum multipole used to infer cosmological parameters. The panels and lines are as in 
Fig. m but pertain to a Stage IV dark energy experiment such as the Large Synoptic Survey Telescope or the European Space 
Agency's Euclid. 



galaxies and galaxy clusters. 

A handful of other recent studies have investigated 
methods for marginalizing over uncertainty in power 
spectra in deri ving co smological constraints from weak 
lensing [l3, [13, [Ull . Refs. [l^,!!!,!!! explore signifi- 
cantly more general parameterizations. However, they all 
reach conclusions that are broadly consistent with ours 
in that each finds self-calibration of uncertainty in the 
nonlinear matter power spectrum a promising approach. 
This broad agreement among different approaches likely 
stems from the well-known fact that cosmological infor- 
mation can be extracted from lensing data based only 
upon geometrical considerations 58, 59^]. Refs. j20l. [57| 
are most similar to ours. These authors explore a halo 
model-based mitigation scheme in which gas and stars 
are modeled separately from dark matter, similar to the 
model proposed by Rudd ct al. 18 1. Semboloni et al. [l^l 



find that their simple model can significantly reduce dark 



energy parameter biases for near term surveys, but that 
improvement may be necessary in order to address Stage 
IV dark energy experiments, a result in broad agreement 
with ours. However, Semboloni et al. [l^ did not use 
their methods to model an independent set of simulations, 
nor did they address the statistical cost of marginalizing 
over additional parameters in their model. 

In summation, our results suggest concentration fit- 
ting as a useful and viable strategy with which to analyze 
cosmological weak lensing power spectra from DES in or- 
der to extract constraints on the dark energy equation of 
state. Based on our analyses. Stage IV experiments may 
remain vulnerable to significant biases in the inferred val- 
ues of the dark energy parameters even after militating 
against baryonic effects with concentration fitting. At 
minimum, estimates for systematic errors such as those 
presented here should be a component of the error bud- 
gets of such experiments. 
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FIG. 8: Biases induced in the dark energy equation of state 
parameter wq from the analysis of a STAGE IV dark energy 
experiment as a function of the maximum multipole used to 
infer cosmological parameters. For STAGE IV experiments, 
four of the OWLS simulations lead to biases significantly 
larger than the others and it is useful to emphasize this. In 
the main panel, the biases that result from analyzing those 
four simulations without any model to account for baryonic 
effects. These biases are clearly very large. Each line is la- 
beled by the name of the corresponding OWLS simulation in 
the panel. The shaded (orange) band, shows the range of bi- 
ases that result after taking the best-fit concentration model 
to describe baryonic effects to analyze these same four simu- 
lations. The biases here are significantly reduced, but remain 
non-negligible (~ Ict). The inset panel shows results for the 
remaining five OWLS simulations. In this case, the inner (or- 
ange) and outer (blue) shaded bands are the same as in Fig. [5] 
for DES. In each of these cases, the mitigation procedure ren- 
ders biases in the dark energy equation of state parameter lOo 
smaller than the statistical error. 



Future work may be able to improve this situation. For 
one, simulations such as the OWLS simulations make pre- 
dictions for the properties of galaxies. It may be possible 
to compare the properties of galaxies in order to deter- 
mine which simulations are more likely to represent the 
observed universe, and use this information to place pri- 
ors on additional parameters in mitigation schemes (the 
concentration parameters in our case, see Ref. (HI). The 
OWLS collaboration has shown that the "AGN" simu- 
lation describes the observed properties of galaxies most 
successfully [13, SI), while our analysis of the "AGN" sim- 
ulation for Stage IV experiments leaves a non-negligible 
residual bias. An important and necessary aspect of fu- 
ture efforts to address these issues with simulations will 
be to develop lensing predictions from baryonic simula- 



FIG. 9: Same as Fig. O but for the bias on Wa from a Stage 
IV experiment. Notice that the vertical ajcis is asymmetric 
about zero. 



tions in larger computational volumes. On another front, 
it may be possible to develop more sophisticated models 
for the influence of baryons on lensing power spectra that 
can minimize biases in inferred cosmological parameters 
without a significant cost in statistical errors. As the cos- 
mological community learns from Stage III experiments 
such as DES and prepares for the Stage IV experiments 
of the coming decade, such efforts should be a high prior- 
ity in order to maximize the scientific yields of the next 
generation of dark energy experiments. 
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